Strong Orientation Effects in Ionization of H2 by Short, Intense, High-Frequency Light Sources 
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We present three dimensional time-dependent calculations of ionization of arbitrarily spatially oriented H J 
by attosecond, intense, high-frequency laser fields. The ionization probability shows a strong dependence on 
both the internuclear distance and the relative orientation between the laser field and the internuclear axis. 
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The ionization dynamics of one- and two-electron pro- 
cesses in diatomic molecules in short, strong laser fields are at 
present under intense experimental investigation (USSl- A 
part of these investigations also focus on the sensitivity of such 
processes to molecular orientation with respect to the light po- 
larization |0]. This is again related to the ultimate goal of con- 
trolling chemical reactions by aligning the reactive molecules 
with respect to each other prior to the intermolecular interac- 
tion n. 

From a theoretical viewpoint such studies are extremely 
complex in the strong-field regime and have been of continu- 
ous interest for nearly two decades (for reviews, see, e.g., [6]). 
In general, only results based on approximate theories such as 
the molecular strong-field approximation |7l [8|] and tunnel- 
ing 0] models have been applied to calculate effects related 
to molecular orientation with respect to the light polarization 
vector. Such approximate theories are, however, often gauge 
dependent IHQiJ] and limited in their applicability to describe 
complex processes. The "slowness" of past and present com- 
puters have, combined with computational challenges related 
to Coulombic multi-center problems, restricted exact theoreti- 
cal calculations including both electronic and nuclear degrees 
of freedom to cases where the internuclear axis is parallel with 
the linear polarization direction Qua or models of reduced 
dimensionality omul. These studies have Riven insight 
into the fascinating interplay between electronic and nuclear 
degrees of freedom; phenomena which at present are beyond 
reach of full dimensional computations. 

In this Letter, we present the first full time-dependent three 
dimensional calculations for the electronic degrees of freedom 
in exposed to a short, strong, attosecond laser pulse. The 
purpose is to follow the behavior of the system with internu- 
clear distance and in particular to display the dependence of 
the dynamics on the angle between the internuclear axis and 
the linear polarization of the field. Calculations are performed 
for 6 cycle pulses with uj = 2 a.u. (23 nm) central frequency. 
This corresponds to pulse durations around 450 as, which 
have already been demonstrated [16]. The ionization prob- 
ability for H(2p) atoms exposed to similar light sources Il7ll 
showed a factor 10 stronger modulation with changing orien- 
tation than what was measured with femtosecond pulses 0]. 
Similar effects in diatomic molecules may thus indicate that 
attosecond pulses may be sensitive probes of the internal nu- 
clear quantum state as well as its orientation. The calculations 



indeed display that the ionization probability depends strongly 
on these parameters. Atomics units (h = m e = e = 1) are 
applied throughout. 

As the nuclear vibrational period is approximately 10 3 
times larger than the pulse duration, the nuclear degrees of 
freedom can be considered frozen during the attosecond pulse. 
Post pulse interplay between nuclear and electronic degrees 
of freedom, which are important for weaker fields, are also 
found to be unimportant here as direct electronic ionization 
dominates. 

The vector potential for the light source is given by 

A(t) = — sin 2 (J t) sm((jt + 0) u p , (1) 

where u p is a unit vector defining the orientation of the lin- 
early polarized field, and <\> is chosen such that the field corre- 
sponding to O represents a physical field [10]. The validity 
of the dipole approximation was investigated in detail very re- 
cently for the present intensity and frequency regime, and was 
found to be well-justified for ionization 11811 . The vector po- 
tential determines the electric field, E(t) = —d t A(t), and the 
translation, ct(t) = J* A(t f ) dt f , which enter the length Hi 
and the Kramers-Henneberger Hkh form fl9ll of the interac- 
tion Hamiltonian, respectively, 
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Hkh = - 



E(*) • r, (2) 



(3) 



|r + R/2 + a(t)| |r-R/2 + a(t)| 



with R the internuclear distance. Both versions of the Hamil- 
tonian have been applied here to secure invariant results. 

For fixed nuclei, we solve the time-dependent Schrodinger 
equation numerically based on a split- step operator approxi- 
mation on a spherical grid. The method was described in de- 
tail elsewhere [SE]]. The wave function is expanded on the 
grid points [(r i: Q jk ) = {r u 6j,<j) k )\ , 

l max ,rnmax 

*(r i ,n jfc ,t;R)= J2 hm{n,t;K)Y Um {Sl jk ) (4) 

l,m 

with origin at the center of the internuclear axis, and with 
parametrical dependence on R. The field-free initial state 



FIG. 1: Ionization probability in the parallel geometry (0 = 0°) 
as a function of the internuclear separation R and the electric field 
strength Eo with to = 2 a.u. and T = 6tt a.u. 



FIG. 2: Ionization probability as a function of the angle between 
the polarization direction and the internuclear axis with uj = 2 a.u., 
Eq — ?> a.u. and T = 6tt a.u. 



|^o) is obtained by the substitution t — > —it in the propaga- 
tor. At internuclear separation R = 2 a.u. this gives an elec- 
tronic ground state energy, e 9 ^ %d = —1.099 a.u., which can be 
compared to the exact value of eo = —1.103 a.u. Reflections 
at the edges r = r max — 60 a.u. are avoided by imposing an 
absorbing boundary. We include up to l max = 15, N r = 1024 
radial points, and in the propagation we use At = 5 x 10 -3 
a.u. Here we stress that for a given lmax, specified at input, 
the quadrature rule of the spherical harmonics uniquely fixes 
the angular points. 

At the end of the pulse, t = T, a fraction of the wave func- 
tion has been removed by the absorber. Since excitation is 
found to be a minor channel, the ionization probability can be 
calculated as P lon = 1 - |(^ |^(T))| 2 . 

The ionization probability as a function of internuclear sep- 
aration and electric field strength is displayed in Fig.[l]for field 
polarization parallel with the internuclear axis. Two striking 
maxima are observed, one for small internuclear separation, 
R ~ 1 a.u., and another for R ~ 3 a.u. When the field 
strength is further increased, the ionization probability de- 
creases, i.e., the molecule is partly stabilized in the intense 
field. This rather counterintuitive mechanism has been studied 
in detail for atoms 1 22] . What happens, is that for increasing 
intensity, the final channels corresponding to single and mul- 
tiphoton ionization close while shake-off dynamics with char- 
acteristic low-energy electron emission becomes the dominat- 
ing ionization mechanism 11811 . Between these two regions a 
"valley" in the ionization curve with increasing intensity may 
occur as indicated around R ~ 1 a.u. and R ~ 3 a.u. for field 
strengths beyond Eo w 7 a.u. At the ground state equilib- 
rium distance, R ~ 2 a.u., and for R ~ 5 a.u., the ionization 
probability is significantly smaller, indicating strong dynamic 
self-interference effects of the electronic charge clouds asso- 
ciated with each scattering center. 

From Fig. [I] we see that the variation in the ionization sig- 
nal is most pronounced for Eo ~ 3 a.u. At this field strength, 
Fig. |2| exposes the ionization probability as a function of in- 
ternuclear separation and as a function of the angle between 



the internuclear axis and the polarization direction of the field. 
An oscillatory behavior of the ionization probability in the 
parallel geometry (0 = 0°) is seen. As increases the os- 
cillations gradually decrease, and in the perpendicular geom- 
etry (0 = 90°), the ionization probability drops monotonically 
with R. In the figure, we also observe opposite functional de- 
pendence with 0: At R ~ 2 a.u. the ionization probability 
increases with while at R ~ 3 a.u. it decreases. 

We now turn to the detailed dynamics and a qualitative un- 
derstanding of the processes underlying the phenomena ob- 
served in Fig. |2| Figure |3 shows snapshots of the wave func- 
tion in the xz-plane at various times for parallel and perpen- 
dicular polarization (the molecule has its internuclear axis di- 
rected along z). In general the photoelectron is ionized in the 
directions of the field. For = 0° the initial charge cloud is 
partly dragged back and forth along the field, and this gives 
rise to a strong interference between various momentum com- 
ponents of the wave function and hence the oscillatory depen- 
dence with R in Fig.[T] This effect is absent at = 90° where 
the two atomic-like charge clouds pertaining to each nucleus 
oscillates in phase back and forth with the electric field. In the 
lower right panel, secondary intensity maxima appear at 30° 
and 150° with respect to the internuclear axis. These struc- 
tures look similar to double slit scattering but have a more 
subtle dynamical origin. 

The oscillatory behavior in Pi on (R) was discussed in a one 
dimensional model earlier in terms of the concept of enhanced 
ionization [23]. Related phenomena were also discussed in 
fast ion-molecule collisions [24]. The following simple ansatz 
gives a qualitative explanation of the oscillations and their ab- 
sence at = 90°: Assume that the outgoing waves from each 
of the scattering centers are a superposition of two outgoing 
spherical waves, 



ife|r+R/2| 

r + R/2| 



ifc|r-R/2| 



R/2| ' 



(5) 



If we take the two scattering amplitudes to be equal, f\ (Oi) = 
JM^). the differential ionization probability can be brought 
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FIG. 3: Tbp: The electric field E(t) of duration T = 6tv a.u. (450 
as) and frequency oo = 2 a.u. The arrows indicate the instants of 
time at which the snapshots of the lower part of the figure are made. 
Color online: Snapshots of the wave function in the xz-plane at times 
corresponding to the beginning (top), the middle (middle) and the 
end (bottom) of the pulse for parallel (left) and perpendicular (right) 
orientation. In all cases the internuclear separation is R = 3 a.u. 
Both the polarization direction and the internuclear axis lie in the 
xz-plane. 



to the form 

dP ion /dn ex 2|/i(fi)| 2 (1 + cos(/cr • R)) 



(6) 



for r ^> R. As seen from Fig. O the main part of the outgo- 
ing wave follows the orientation of the field. Hence we expect 
that for parallel polarization, the main contributions will be 
for f parallel to R, which again gives raise to oscillations in 
R with wave number k. Given that the one-photon ionization 

dominates, we find the wave number as k 



f 2(uj-I p (R)), 

where I p is the effective ionization potential. This is seen to 
be consistent with the results in Figs, [l] and and we have 
also confirmed these findings for other values of uj. The ab- 
sence of oscillations in the case for perpendicular polarization 
is understood accordingly: The wave is sent out mainly in the 
direction given by = 90 °, and since the outgoing waves will 
have no phase difference due to the separation of the scatter- 
ing centers in this direction, no interference pattern will occur 
(f R = 0). 



FIG. 4: Angular photoelectron spectrum in the scattering plane for 
parallel (dashed curve) and perpendicular (full curve) geometry as 
a function of the polar angle for a 6 cycle field with Eo = 3 a.u. 
and uj — 2 a.u. In both cases, the angle denotes the direction of the 
outgoing electron with respect to the internuclear axis. We have used 
equal normalizations for the two curves. 



The monotonous decrease in Pion with R at 6 — 90° is 
reasonable since the decrease in the ionization potential leads 
to an effective higher final state electronic momentum with 
increasing R. 

The angular distribution of the ionization probability can be 
calculated from the time integral of the radial current density 
through the solid angle element dVt at a chosen distance a 
from the origin 



dRc 



dn 



dt j(a, t) • f 



dtSm ( — 
or 



(7) 



where the distance a is chosen large enough to exclude con- 
tribution to the current from the quiver motion of an electron 
close to the nucleus and small enough to avoid effects induced 
by the absorber. The application of this procedure to the out- 
going waves of the lower panel of Fig. [3] results in the inten- 
sity spectrum of Fig.|4| As already observed from Fig. [3] the 
photoelectrons are most likely ejected in the direction of laser 
polarization. In the perpendicular arrangement the two peaks 
around 30° and 150° can also be quantified: The integrated 
probability connected to these secondary maxima amounts to 
about 25% of the total ionization probability. As these peaks 
originate from interference of the outgoing waves, it is inter- 
esting to note that in a very recent calculation of High Har- 
monic Generation (HHG) in a reduced model with respect to 
the electronic degrees of freedom such interference does not 
occur l25ll . In that work it is pointed out that orientational 
effects are very important for HHG. 

In conclusion, fully non-perturbational calculations of the 
ionization dynamics of H^ molecules in intense attosecond 
light sources have been carried out. Very strong orientation 
effects have been found, demonstrating that in order to obtain 
a full understanding of the molecular ionization dynamics, all 
three electronic degrees of freedom must be included. The 
qualitative features are determined by interference effects re- 
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lated to double-center scattering and the distinct features in 
the electron spectra show that intense attosecond pulses can 
resolve the instantaneous vibrational and orientational quan- 
tum state of diatomic molecules. 
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